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Abstract 

In this paper, we are interested in approximating the solution to scalar conservation laws using 
systems of interacting stochastic particles. The scalar conservation law may involve a fractional 
Laplacian term of order a £ (0,2]. When a < 1 as well as in the absence of this term (inviscid 
case), its solution is characterized by entropic inequalities. The probabilistic interpretation of the 
scalar conservation is based on a stochastic differential equation driven by an a-stable process 
and involving a drift nonlinear in the sense of McKean. The particle system is constructed by 
discretizing this equation in time by the Euler scheme and replacing the nonlinearity by interaction. 
Each particle carries a signed weight depending on its initial position. At each discretization time 
we kill the couples of particles with opposite weights and positions closer than a threshold since the 
contribution of the crossings of such particles has the wrong sign in the derivation of the entropic 
inequalities. We prove convergence of the particle approximation to the solution of the conservation 
law as the number of particles tends to oo whereas the discretization step, the killing threshold 
and, in the inviscid case, the coefficient multiplying the stable increments tend to in some precise 
asymptotics depending on whether a is larger than the critical level 1. 

Introduction 

We are interested in providing a numerical probabilistic scheme for the fractional scalar conservation 
law of order a 

d t v(t,x) + a a (-A)%v(t,x) + d x A(v(t,x)) = 0, (i,x)eR+xR, (0.1) 

where — (—A) "2 is the fractional Laplacian operator of order < a < 2 (defined in Section 2), and A 
is a function of class C 1 from R to R. We also consider the equation obtained by letting a — > 
in (|0.2|) . namely the inviscid conservation law 

d t v(t,x) + d x A(v(t,x)) = 0, (i,i)el + xl. (0.2) 

In [9] [10], these equations are interpreted as Fokker-Planck equations associated to some stochas- 
tic differential equations nonlinear in the sense of McKean, which can be approximated by a particle 
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system. We introduce an Euler time discretization of this particle system and show the conver- 
gence of its empirical cumulative distribution function to the solution of (|0.ip . We also study its 
convergence to the solution of (|0.2p as the parameter a goes to 0. 

Euler schemes for viscous conservation laws have already been studied in [3], [4], [5] or [6], where 
a convergence rate of -7= + \TKt is derived in the case a = 2, N denoting the number of particles, 
and At being the time step. 

To give the probabilistic interpretation to (|0.1[) we consider the space derivative u = d x v of a 
solution v to equation (|0.1I) . which formally satisfies 

dtu t = —a a (—A)%u t - 8 X (A'(H * u t )u t ) , (0.3) 

where H — l[o,oo) denotes the Heavyside function. When uo is a probability measure, that is, when 
the initial condition vo of Equation (|0.ip is a cumulative distribution function, Equation (|U.3f) is the 
Fokker-Planck equation associated to the following nonlinear stochastic differential equation 

(dX t =adL? + A'(H*ut(X t ))dt 
\u t — law of X t 

where Lf is a Markov process with generator — ( — A) 2", namely v2 times a Brownian motion for 
01 — 2, and a stable Levy process with index a in the case a < 2, that is to say a pure jump Levy 
process whose Levy measure is given by c a r^Wfs ; where c a is some positive constant. 

We can still give a probabilistic interpretation to Equation (|0.1[) if the initial condition vo has 
bounded variation, is right continuous and not constant. Indeed, in that case vo can be written 
as vo(x) — a + j^^duoiy) — a + H * uo(x) for some finite measure wo- By replacing vo(x) by 
(v (x) - a) Ool(R))- 1 and A(x) by A(a + a;|uo|(R))(|iio|(R)) _1 in 63) (K denoting the total 
variation of the measure uq,), one can assume without loss of generality that a — and that |uo| is 
a probability measure. We denote by 7 = jfj^j the Radon-Nikodym density of uo with respect to 
its total variation. Notice that 7 takes values in {±1}. 

Then, Equation (|0.3p is the Fokker-Planck equation associated to 

(dX t =adL? + A'(H*Pt(X t ))dt 
\P = law of X 

where P denotes the measure defined on the Skorokhod space D of cadlag functions from [0, 00) to R 
by its Radon-Nikodym density ^-p — j(f(0)), with / the canonical process on V, and Pt denotes its 
time marginal at time t, i.e the measure defined by Pt(B) — 7(/(0))ls(/(t))dP(/), for any B 
in the Borel a— field of R. 

The rest of the paper is organized as follows: 
In Section 1 we define the particle approximation for the stochastic differential equation (|0.4p . 
Section 2 is devoted to the definition of the different notions of solutions used in the article. 
In Section 3, we analyze the convergence of the time-discretized particle system to the solution of 
the conservation law in different settings : for both a constant or vanishing diffusion coefficient and 
any value of < a < 2. 

Finally, we present some numerical simulations in Section 4. Those simulations are compared with 
the results of a deterministic method described in [7] . 

In the following, the letter K denotes some positive constant whose value can change from line 
to line. 

Acknowledgement: We warmly thank Jerome Droniou for providing us with the article [7] 
and the source code from the corresponding deterministic numerical scheme, and Mireille Bossy for 
fruitful discussions about the convergence rate of the numerical schemes. 
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1 The particle approximation 



In this section we construct a discretization of (|0.4[) consisting of both a particle approximation 
in order to approximate the law of the solution and an Euler discretization to make the particles 
evolve in time. The idea is to introduce N particles X N '* , . . . , X ' which are N interacting copies 
of the stochastic differential equation (|0.4p . where the actual law P of the process is replaced by the 
empirical distribution of the particles i> J^ 4=1 S x iv,i. 

In continuous time, those particles are driven by JV independent Brownian motions or stable Levy 
processes with index a and undergo a drift given by A' (H * j$ (.)) , with ft^ = i J3 j=1 r y(X^ r ' x )S x N,i . 

The natural way to introduce the measure fit in the dynamics is to give each particle a signed weight 
equal to the evaluation of 7 at the initial position of the particle. Then, H * jlf (x) is simply given 
by the sum of weights of particles situated left from x. 

The entropy solution to (10. 1[) has a non-increasing total variation (see [2]), which can be inter- 
preted probabilistically as a compensation of merging sample paths having opposite signs. For a 
more precise statement in the case a — 2, see Lemma 2.1 in [9]. It is thus natural to adapt this 
behavior in our particle approximation by killing any merging couple of particles with opposite signs. 

In [9] Jourdain proves, for a = 2 in continuous time, the convergence of the particle system to 
the solution of the nonlinear stochastic differential equation through a propagation-of-chaos result. 
Moreover, the convergence of the signed cumulative distribution function H * p,^ to the solution 
to Equation (|0.1|l is also proved, as well as convergence to the solution to the inviscid equation as 
(7 — ^ 0. In [lOj the same results are generalized to the case 1 < a < 2, assuming 7 = 1 in the case 
of a vanishig viscosity. However, to our knowledge there is even no existence result for the particle 
system in continuous time when a < 1, since the driving Levy process is somehow weaker than the 
nonlinear drift. 

In discrete time, the probability of seeing two particles actually merging is 0. To adapt the 
murders from the continuous time setting, we thus kill, at each time step, any couple of particles 
with opposite signs separated by a distance smaller than a given threshold ejv going to zero as N 
goes to 00. Though, one has to be careful, since one can have more than two particles lying in a 
small interval of length en- Precisely, the particles are killed in the following way: kill the leftmost 
couple of particles at consecutive positions separated by a distance smaller than the threshold £jv 
and with opposite signs. Then, recursively apply the same algorithm to the remaining particles. 
This can be done with a computational cost of order O(N). The essential properties satisfied by 
this killing procedure are the following: 

• to each killed particle is attached another killed particle, which has opposite signs and lies at 
a distance at most ejv of the first particle. 

• after the killing there is no couple of particles with opposite signs in a distance smaller than ejv. 

• the exchangeability of the particles is preserved. 

• after the murder, the quantity H * (X^' 1 ) remains the same for any surviving particle. 

We are going to describe the killed processes by a couple (/, «) in the space K, = T> x [0, 00] of 
cadlag functions / from [0, 00) to R endowed with a death time k £ [0, 00]. The space K, is endowed 
with the product metric d((f,Kf), (g,K g )) = ds{f,g) + | arctan(fc/) — arctan(ft 9 )| , where ds is the 
Skorokhod metric on D, so that (IC,d) is a complete metric space. It could seem more natural to 
consider the space X>([0, 00), RU {d}) of paths taking values in R endowed with a cemetery point d. 
However the corresponding topology is too strong to prove Proposition 13.41 

The precise description of the process is the following: each particle will be represented by a 
couple {X N '\^) G JC. Let (X^) i£N be a sequence of independent random variables with common 
distribution |wo| and let /ijv > denote the time step of the Euler scheme. At time 0, kill the 
particles according to the preceding rules, that is to say, set nf = for killed particles, which will 
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not move anymore. Those particles will not be taken into account anymore. Now, by induction, 
suppose that the particle system has been defined up to time fc/ijv, and kill the particles according 
to the preceding rules (i.e. set nf = fcftjv and X^'' = X kf f for all t > fc/ijv, if the particle with 
index i is one of those). Then let the particles still alive evolve up to time (fc + l)/ijv according to 



dXf' 1 = A' 



N z * ^kh /v —^khj\ 



where (Z/); s n is a sequence of independent a-stable Levy processes for a < 2, or a sequence of 
independent copies of y/2 times Brownian motion, which are independent of the sequence (Xq)^^. 
The particle system is thus well-defined, by induction. 

Let ji N = Si=i ^{x N - i k n ) ^ be the empirical distribution of the particles. For a 

probability measure Q on K, and t > 0, we define a signed measure Qt on R by: 

Qt(S)= / l s (/(t))l K>t7 (/(0))dQ(/ J «), 

JtC 

for any B in the Borel cr— field of R. With these notations, on the interval [fc/tjv, (fc + l)/ijv), a 
particle, provided it is still alive, satisfies 



AX 



N,i 



A' (H*H^ S (x^))dt + a N dL\. 



Notice that the sum of the weights of alive particles fi^ 1 (R) = ^2 K N >t j(Xq) is constant in time, 
since the particles are killed by couples of opposite signs. 



2 Notion of solutions 

In this section, we recall the different notions of solutions that are associated to the equations (|0.1[) 
and (|0.2[) . Indeed, due to the shock-creating term d x (A(u t )), the notion of weak solution is too 
weak, and does not provide uniqueness when the diffusion term is not regularizing enough. The 
best suited notion in those cases is the notion of entropy solution. 

In [11], Kruzhkov shows that for vo £ L°°((0,oo)) existence and uniqueness hold for entropy 
solutions to H0.2|) . defined as functions v £ L°°((0, oo) x R) satisfying, for any smooth convex 
function 77, any nonnegative smooth function g with compact support on [0, 00) x R and any ip 
satisfying tp' = rj'A', the entropic inequality 

r){v )g + J (J n(vt)d t gt + il){v t )d x gt\ dt > 0. (2.1) 

It is well known that this entropy solution can be obtained as the limit of weak solutions to (|0 . 1 f) 
as a — > in the case 01 — 2. 

Weak solutions to (|0.1|) (see [9]) are defined as functions v € L°°((0, 00) x R) satisfying, for all 
smooth functions g with compact support in [0, 00) x R, 

roc p f" 00 p p 00 f 

^o<?o + / / v t dtg t dt-a a / / v t (-A)%g t dt+ / / A(v t )d x g t dt = 0. (2.2) 
Jo in Jo Jr Jo Jm 

For a < 2, we denote by (—A) 2" the fractional symmetric differential operator of order a, that 
can be defined through the Fourier transform: 
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An equivalent definition for (—A) 2 uses an integral representation 



(-A) 2 u(x) 



u(x + y) - u(x) - l M < r u'(x)y 

\y\ 1+a 



dy 



for any r 6 (0, 00) and some fixed constant c a (see [8]), depending on the definition of the Fourier 
transform. 

It has been proven in [9] and [10] that existence and uniqueness holds for weak solutions of (10.11) . 
for 1 < a < 2. However, for < a < 1, the diffusive term of order a in (|0.ip is somehow 
dominated by the shock-creating term, which is of order 1, so that a weak formulation does not 
ensure uniqueness for the solution. We thus have to strengthen the notion of solution, and use 
entropy solutions to (|0.ip . defined in [5] as functions v in L°°((0, 00) x R) satisfying the relation 



for any r > 0, any nonnegative smooth function g with compact support in [0, 00) x R, any smooth 
convex function r\ : R — > R and any ip satisfying ip' = rj'A'. Notice that from the convexity of 77, 
the entropic formulation (|2.3|) for a parameter r implies the entropic formulation with parameter 
r' > r. Also notice, using the functions T](x) = ix that an entropy solution to ()0.1[> is a weak 
solution to (|0.ip . 

In [2], Alibaud shows that existence and uniqueness hold for entropy solutions of (JOTTJ pro- 
vided that the initial condition vo lies in L°°(R). The entropy solution then lies in the space 
C([0, 00), L 1 ( jj^s )). He also proves that the entropy solution to (|0.ip converges to the entropy 
solution to (HEU) in the space C([0, T], Lj^R)) as a -> 0. 

3 Statement of the results 

The aim of this article is to prove the three following convergence result, each one corresponding to 
a particular setting. 

Theorem 3.1. Assume < a < 1. Let crjv = a be a constant sequence. Let en and Iin be two 

vanishing sequences satisfying the inequalities 




(2.3) 



N x < 4 sup \A'\h,N < £n, and N 1 
[-1,1] 




for some positive X. For a = 1, ako assume /ijv < ejvA' 



-i/A 



. A holds for any T > 0, 




where vt denotes the entropy solution to the fractional conservation law (|0,1[) • 
Theorem 3.2. Let en, hM and o~n be three vanishing sequences such that 



sup |^4.'j/ijv < en 
[-1,1] 
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for some A > 0. If a > 1, also assume on <e„ °W A . Then, for any T > 0, 

fT .. .. 
lim / Epftjif-D, . N d£ = 0, 

N^ovJ II ,lLl (rftO 

where Vt denotes the entropy solution to the inviscid conservation law (|0.2|l . 

The additional assumption for a > 1 comes from the fact that in this case, the dominant term 
is the diffusion, while in the limit there is no diffusion anymore. The assumption ensures that the 
diffusion is weak enough not to perturb the approximation. For a < 1, the dominant term is the 
drift, as in the limit, so that no additional condition is needed. 

Theorem 3.3. Assume 1 < a < 2. Let <tjv = a be a constant sequence, and let en and fejv be two 
vanishing sequences. It holds for any T > 0. 

rT ., 

lim / E\\H* ur -v t \\ , , dt = 0, 

jv^ooy ii pt '^(t^t) 

where Vt denotes the weak solution to the fractional conservation law (10. 1[) . 

In order to prove those three theorems, we will have to control the probability of seeing particles 
merging. In the case a < 2, this is mainly due to the conjunction of the small jumps of the stable 
process and the drift coefficient, while the large jumps of the stable term do not play an essential 
role. As a consequence, for a < 2, we consider another family of evolutions coinciding with the 
Euler scheme on the time discretization grid, for which we consider differently the jumps which are 
smaller or larger than a given threshold r. The choice of this parameter has to be linked to the 
parameter r appearing in the entropic formulation (|2.3|l . since they play a similar role: the third 
term in (|2,3|l corresponds to the effect of jumps larger than r in the driving Levy process and the 
fourth term corresponds to jumps smaller than r. This evolution is designed so that on the first 
half of each time step, the process will evolve according to the drift and the small jumps, and on 
the second half of each time step, it will evolve according to the large jumps. More precisely, let 

v\&y,At) = J2 <Wj,t) 

be the jump measure associated to the Levy process L 1 and let 

dydt 



v\Ay,At) = v l (dy,At) - c a - 
be the corresponding compensated measure, so that 



L\ = I yv\Ay,At)+ I yv{Ay,At), 

J(0,tlx-flul>r> J fO.«lx-flul<r> 



'(0,t]X{|s/|>r} -I (0,t}x{\y\<r} 

where the right hand side does not depend on r. We define the process X N ' l ' r by 

x N,i,r = xl + o N L N ' % ' r + o N A N - l ' r + A N '\ 

where 



• L t ' l ' r is the large jumps part defined by 



yi>*(Ay,As), 

(0,o(t)]x{| H |>r} 
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, MX \kh N for te [kh N ,(k + l/2)h N ] 

where alt) = < . Ihis process is 

[kh N + 2(t-(k + l/2)h N ) forte [(k + l/2)h N , (k + l)h N ] 

constant on intervals [khN, (k + l/2)/ijv] and behaves like a Levy process with jump measure 

l| B |> r ofri on intervals [(fc + l/2)ftjv, (fc + l)/ijv]. 

• Af r,1 ' r is the small jumps part, defined by 

Af' l ' r = / £ l (dy,ds), 

i(0,i,(t)]x{|i/|<r} 

( kh N + 2(t - kh N ) ior te[kh N ,(k+l/2)h N ] 
where bit) — < . Ihis term behaves like 

\(k + l)h N for te[(k + l/2)h N ,(k + l)h N ] 

a Levy process with jump measure l| y |<r ofji on intervals [khN, (fc + 1/2)/ijv] and is constant 

on intervals [(k + 1/2)Hn, (k + l)/ijy]. Notice that the process A^' 1 ' 1 " is a martingale. 

• yl^' 8 is the drift part, which satisfies A^'' = 0, is constant over each interval [(fc + l/2~)hu, (k + 
l)/iiv], and evolves as a piecewise affine process with derivative 2A'(H * fikh N (^khlr)) on 
intervals [khN, (k + l/2)/ijv]- 

One can check that for any r, the process (X N ' 1 ' r , . . . , X N ' N ' r ) is equal to (X N '\ . . . , X N ' N ) 
on the time discretization grid up to killing time. Conditionally on the positions of the parti- 
cles at time khN, the particles evolve independently on [khN, (k + l)/ijv], and the evolution on 
[khN, (k + 1/2)/ijv] is independent of the evolution on [(k + 1/2)/ijv, (k + 1)/ijv]. Since the entropic 
formulation (|2.3[) with parameter r is stronger than the one with parameter r' > r, we have to make 
the parameter r tend to zero in order to prove the entropic formulation for any parameter. However, 
this convergence has to satisfy some conditions with respect to N, ft, at and ejv. We will explain later 
why a suitable sequence rjv exists under the conditions given in the statement of Theorem 13.11 

In order to prove Theorems 13.11 and 13.21 we introduce fi N ' r the empirical distribution of the 
processes {X N ' x,r 



and by 7r ' r the law of /x' 



t N,r 

The following proposition is the first step in the proof of Theorems 13.11 13.21 and 13.31 

Proposition 3.4. • Assume a < 2. For any bounded sequences (/ijv), (cjv) and (ejv), and for 
any sequence (rjv), the family of probability measures (tt ' rjv )jveN is tight in V(V(fC)). 

• Denote by ir N the law of fi N . For any bounded sequences (/ijv), (cjv) and (ejv), the family of 
probability measures (7r JV )jvgM is tight in V (V (K.)) . 

Proof. We first check the tightness of the family (ir N ' TN )jvsn- 

As stated in |13| . checking the tightness of the sequence ir N > rN boils down to checking the 
tightness of the sequence (Law(X JV,1,rjv , k^)). Owing to the product-space structure, we can check 
tightness for X N ' 1,rN and Ki separately. 

Of course, tightness for /tf is straightforward since it lies on the compact space [0, oo], and it is 
enough to check tightness for the laws of the path (X N ' 1,rN ). For simplicity, we will assume that 
A — 0, which is not restrictive since A' is a bounded function so that the perturbation induced 
by A belongs to a compact subset of the space of continuous functions, from Ascoli's theorem (also 
notice that the addition functional from PxC([0, oo)) to T> is continuous). We use Aldous' criterion 
to prove tightness (see [T]). First, the sequences (X^' 1 '' N )jvsn and (sup[ T i | AX^' 1 ' 1 "™ |)jvgN are 
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tight, since (X^' 1 ^™) is constant in law and ^sup[ T j j AX N,1 ' rN | j is dominated by the identically 
distributed sequence 

AL 1 



Sup (J at I SUp 

N } [0,T+su Pjv h N ] 



N 



v N,l,r N v N,l,r N 



where 



Then let tn be a stopping time of the natural filtration of X ' ,rjv taking finitely many values, 
and let (5n)n€N be a sequence of positive numbers going to as TV — > oo. One can write 

>s) <P(^|A^- -A^|> £ /2) 

+ r(a N \ L* 1 ™ - L?f» | > e/2) (3.1) 

sup o N \Cf rN \ > e/2] + P ( sup aAr|£ t >rN | > e/2 ] , 
ts[o,s N ] J \te[o,s N ] J 

Cf r = j yv{dy,dt) and Cf r = / yu(dy,dt), 

J (0,t]x{\y\<r} J (0,t]x{\y\>r} 

the measure v being the jump measure of some Levy process C with Levy measure ^1+1 , and & is 
the compensated measure of v. Now, using the maximal inequality for the martingale (£f rN )te[o,s N ] > 
noticing that (£.g'^)re[o,i] m ^ so a martingale, we deduce 

sup \Cf rN \ > e/2o N ] < sup P[ sup \Cf r \>e/2o N 

ie[0,<5 N ] J r6[0,su Pjv r N ] \t£[0,S N ] 

<2on£ 1 sup E([£j-' 

r£[0,sup N rjv] 

= 2a Jve - 1 E(|/:f; uPNrN | 

— > 0. 

JV— >oo 



For the large jumps parts, one writes, 

sup \Cf rN \ > e/2a N ) < P f sup \C t \ + sup \Cf rN \ > e/2a N 

sup \Ct\ > e/4criv ] +P [ sup \Cf rN \ > e/4a N 
te[o,s N ] j yteio.a^r] 

-> 0. 

iV— >co 

As a consequence, the family (Law(X JV,1 ' rjv ))jvgN is tight in T>. 
Thus, the family (7r JV ' rN ) JveN is tight. 

The proof is essentially the same for the tightness of (7r JV )jvgN, with a few simplifications, since 
we do not treat separately large and small jumps. It also adapts in the case a = 2, since the 
Gaussian distribution has thinner tails than the a— stable distribution for a < 2. □ 

The use of the path space K, instead of 2?([0, oo),RU {d}) for a cemetery point d is crucial in 
the proof of Proposition 13.41 since in the latter case, we need to control the jumps occuring close to 
the death time in order to prove tightness. The following example is illustrative: if we consider a 
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sequence /„ of paths starting at 0, jumping to 1 at time 1 — 1/n, and being killed at time 1, then /„ 
does not converge in 2?([0, oo),R U {d}), while it does in K,. 

The following lemma deals with the initial condition of the particle system. 
Lemma 3.5. If it 00 is the limit of some subsequence of ir N or ir N ' rN , then for 7r°° -almost all Q, 
for all A in the Borel a— field o/R, 



Qo(A):= / l K>0 l meA dQ(f,n) = \uo\(A). (3.2) 
Jm 

In particular, k is Q— almost surely positive for 7r°°-almost all Q. 

Proof. In a first time, we control the probability of seeing a partincle dying within a short time. 

Let us write the Hahn decomposition Uq — of the measure uo, the measures Uq and Uq 
being positive measures supported by two disjoint sets B + and B~ . From the inner regularity of the 
measure Ug , for any 8 > 0, one can find a closed set F + C B + such that Uq(F + ) > Ug(B + ) — 8. 
The complement set 0~ = (F + ) c is then an open subset of R, which can thus be decomposed as a 
countable union of disjoint open intervals 0~ = Um=i] am '^ m [- For a large enough integer M, and 
for es > small enough, the set O 6 = U^Jnm + £s,b m — eg[ is such that Uq(O s ) > Uq(0~) — 8. 
Consequently, we can write R as a partition 

R = F + U (B~~ n O s ) U B s , 

where B s = (F + U (B~ n 6 )) c has small measure |ito|(0*) < 28, particles starting in F + have a 
positive sign, and particles starting in (B~ n O s ) have a negative sign. Let N be large enough to 
ensure ejv < Eg /3. The distance between any element of F + and any element of O 5 is larger than eg. 
As a consequence, if the particles with index i and j kill each other before time r, then either one 
of them started in B s , or one of the particles i and j moved by a distance larger than £5/3. This 
writes 

|J {i, ref < Tg} = 2J < j, X N ' 1 '™ and X N ' J ' rN kill each other j 



< 2(t ^ i,X™' l ' r ~ G 23° or sup \X^ l ' rN - X?,\ > eg/3 

{ t£[0,r] 

As a consequence, if ts > is small enough so that P(sup te[0:Ti5] \x"' l ' rN - X£ '*' rjv | > eg/3) < <5 
(this can be achieved using an adaptation of (|3.1[l ). it holds 

< **) = (j [i, «f < rg}) < h) |i, S B s or jup ^ |xf Ap " - Xj| > e^J j 

< 2P (xj G B s ) + 2P f sup |Xf ■ < ' rjv - Xf ' i,rjv I > eg/3 

v y \te[o,r 6 ] 

< 68. 

Consequently, 

E 7r °°(Q(K < t)) < liminf E"'™ (Q(k < r)) = liminf P(«f < r) -> 0. 

JV JV T->0 

Thus for 7r°°-almost all Q, k is Q— almost surely positive. As a consequence, for any bounded 
continuous function tp, 



l«>o*»(/(0))dQ(/c,/)- / ^d| Mo 



E^ 



/ *»(/(0))dQ(K,/)- / pd|uo|| 

^"1/ ^(/(0))dQ( K ,/)- / pd|, 
\Jk Jk 
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from the law of large numbers. 

□ 



The main step in the proof of Theorems 13.11 [3721 and 13731 is the following proposition: 
Proposition 3.6. Let en and /ijv be vanishing sequences. 

• If ajv is a constant sequence and < a < 1, suppose iV~ 1/,Q < N~ 1 ^ x en and N~ x < 
4supj_ 1 jj \A'\1in < en for some positive A. If a = 1, also assume ftjv < N~ 1 ^ x en- Then, 
there exists a sequence (rjv) of positive real numbers, such that the limit of any converging 
subsequence of -k n ' Tn gives full measure to the set 

{Q G V(JC),H * Qt(x) is the entropy solution to (|0.ip }. 

• Let on be a vanishing sequence and assume N~ x < 4supr_ 1 \A \Hn ^ en for some positive A. 
If 1 < a < 2, also assume crjv < £ N a iV~ 1//A . Then 

{Q G V(JC),H * Qt(x) is the entropy solution to (|0.2|l I 

is given full measure by any limit of a converging subsequence of tt ' rN , for a well chosen 
sequence (rjsr), in the case a < 2, and by any limit of a converging subsequence of tv n if a — 2. 

• If crjv is a constant sequence and 1 < a < 2, t/ie Zimit o/ any converging subsequence of tt n 
gives full measure to the set 

{Q G V{K),H * Qt(x) is the weak solution to (10.1[) }. 



Proposition 13.61 will be proved in Section 13.11 We first admit it to end the proofs of Theo- 
rems 13.11 13.21 and 13.31 We need the following lemma. 

Lemma 3.7. Let a < 2 and rN be a sequence of positive numbers going to zero. Then it holds, for 
any T > 0, 

Jim / E\\H*fi,? -H*ji?> r »\\ , dx xdi = 0. 
Proof. It holds, by exchangeability of the particles, 



K, ■ >t 

<j\(l K ^ t \x^-Xf^\A^) 



I d:rdt 
<*l x 2 + 1 

dt. 



This last quantity goes to zero, since the processes X 1 ^' 1 and X iY ' 1 ' rjv coincide on the discretization 
grid, whose mesh vanishes. Indeed, for t G [fc/ijv, (k + l)/ijv) 

E (l K „ >t \X?^» _ X AU, A7r ) < E ( 1<>t | X -.i^ _ | A7r ) +E >t |Xf - 1 - Xf4»| A tt) 

<Kh]( 2 . (3.3) 
For this last estimate, we used, for an a— stable Levy process L, the inequality 

E(|L t |Al)<iflE(j.L t | a!/2 ) =Kt 1/2 . 

□ 
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From Lemma 13771 it is sufficient to show lim j' K\\H * 

JV^oo J0 " 1 

Theorems ED and [5~2l 



ftlll^i^ da. )dt in order to prove 



Proo/ 0/ Theorems \3.1H3. e 3373\ We write the proof for Theorems 13 . 1 1 and 13721 in the case a < 2. The 
proof of Theorem 13.21 with a — 2 and Theorem 13.31 is the same, with 71-^ replacing 7r JV ' rjv . 

Let j k be a Lipschitz continuous approximations of 7, with P(7(Xq) 7^ 7 fe (Xo)) < 1/fc (see [9], 
Lemma 2.5, for a construction of such a 7*). We have, by exchangeability of the particles, 



E 
<E 
+ E* 



Fl. 

JO i/M 



H*ij^' r "{x)-vt(x)\ 



dxdt 



kr v N,l,r* 



dxdt 



l K>t H(x - f(t))y K (/(0))dQ(/, k) - v t (x) 



X 2 + 1 

da; 



(3.4) 



x 2 + 1 



From the assumption on 7 fe , the first term in the right hand side of (|3.4[) is smaller than 2Tv/k which 
vanishes as k goes to 00. The bounded function 



l K>t H(x - /(t))7 fe (/(0))dQ(/, k) - vt(x) 



da;dt 
x 2 + 1 



is continuous. From Proposition 13.61 the second term in the right hand side of (|3.4[) converges, as N 
goes to 00 to 



l K>t H(x - f(t)) 7 fc (/(0)) - 7(/(0)) dQ(f, k) 



Ax 



x 2 + l 



This terms goes to zero as k tend to infinity using the argument of the begining of the proof 
with X N,1 ' rN replaced by the canonical process y. 

□ 



3.1 Proof of Proposition [3761 



This section is devoted to the proof of Proposition 13.61 Since the hardest part of this proof is the 
first two items, we do not give all details for the third item and for the second one in the case a = 2. 
Indeed, for these two last settings, the separation of small jumps and large jump is not necessary 
for the proof. 

Let rjv be a sequence of positive real numbers, going to zero as N — > 00, which will be explicited 
later. Let r > and c be reals numbers, r\ a smooth convex function, ip a primitive of A'rf and g 
a smooth compactly supported nonnegative function. We define the function <pt{x) = f_ gt{y)dy. 
Note that <p is smooth, and nondecreasing with respect to the space variable. We consider a 
subsequence of w ' rN , still denoted ir N ' rN for simplicity, which converges to a limit ty°° . We want 
to prove that, for 7r°°— almost all Q, the function H *Q t satisfies the entropy formulation associated 
to the corresponding case. 

One can write, for any k > and t £]fc/ijv, (k + l)/ijv] 



/=]•„■ N . N _ , v N,i,r 

I K t A Kj > t,X t 
\ I P I 3i,j, A Hif > t, ijnZ^ 



E 



-1 • ■ N . N . , v N,i,r 

Hi A Kj > t,X t 



i,j,k,N 



= x 



X N ' 1 + A N ' j - A N ' 1 

A kh N ' ■ / H 



(Xkhl,)l 
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where we denote 



i,j,N,k j^N,i,r ^N,i,r ^N,j,r j^N,j,r j_ j N,i.r j N,i,r jN,j,r j_ jN,j,r 



A 



kh \ 



J kh a 



+ L 



k h ^ 



From the conditional independence of the processes L N ''' r , L N '-'' r , A N ''' r and A JV '-'' r , the random 
variable Z\'^' ' has a density. As a consequence, since the process A^^ — A 1 ^' 1 is deterministic 
on [kliN, (k + l)/i]v] conditionally to (A^'^)^, the above probability is zero, meaning that for all 
time t > 0, the alive particles X t JV,I ' r,v almost surely have distinct positions. As a consequence, the 
function r\ (H * jl t ' rN (x) ) is the cumulative distribution function of the signed measure 



£ = E ^ 



where 



V 



i 

N 



E 



7^) 



V 



1 

N 



E y&i) 



JV<x JV,,,r JV 



= i Kf >t („ (h * Af' rN (xf — )) - r, (h * p? 

Let (Cm)m6N be the increasing sequence of times which are either a jump time for some L JV,l,rjv (i.e. 
a jump of size > tn for X N,l ' rN ) or either a time of the form fc/ijv/2. One has 

oc 

-(to,<po) = (tcLi'Pc™) - (^cl-l^Cm-l) 

771= 1 

OO 

= E E (*>C» (Af m ^) - (xf TO --)) (3.5) 

K iV >0 m = l 
oo 

+ E E K^, «r*) -wi m _ ltp(m (xg 



ref >0 m =! 

Notice that these infinite sums are actually finite, since the function ipt is identically zero when t is 
large enough, and since the process 

( L N,i,r N ^ l n,n,t n -j has a finitg number of 

jumps on bounded 

intervals. 

We consider the first term in the right hand side of (|3.5p . Denote by v % ' r = < ^(Ai, Jv i '*"+AA JV 

the jump measure associated to L N ' z ' r + A > t ' r , and by 



v*' r (dy, dt) = v l ' r {dy, dt) - 2c a {x*^\y\<r + (1 - xf)l|y|>r) j^ dt 



its compensated measure, where Xt 1 = Xl^Lo l[fch N ,(fc+i/2)h N ) W- Let us apply Ito's Formula on the 
interval (Cm-iiCm). If Cm-i = &/Jiv for some integer fc, then £ m = (k + l/2)/ijv, and almost surely 
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= X^'*'t., holds. As a consequence 

(k+±)h N - (fc+i)h N 

/ v N,i,r \ ( v N,i,r\ 

(k+^)h N f(k+^)h N 

d m {X^)At + 2 / d.vt{X?> i <*)A! (H * (Xf^)) At. 



/ / [<Pt{Xt_l %lT + a N y) - ip t (X^' r ) - a N yd x ifit(X^ l ' r )) v l ' r {Ay,Ai) 

J (kh N ,(k+l/2)h N ) J {\y\<r} V 

/ d^ix^l l f yv' r (dy,dt) 

J (kh N ,(k + l/2)h N ) \J{\y\<r} J 



+ 

J(kh N ,(k + l/2)h N ) \J{\y\<r} 

If (m-i is not of the form fc/ijv, then the process X N ' l ' r is constant on the interval [Cm-i,Cm), 
and one has tpQ m {X^f) — LpQ m l (X^ r '^ i ) — /5 m dt'fit{X^ I ' l ' r )At. Summing over all the intervals 
(Cm-ijCm), Equation (|3.5|) writes, denoting r t = max{( m ,C m < t}, 

-(^o) = E f°< (9 t ^(Xf^«) + 2 x fa^^Xf^")A'(i/*An' rN (^' rN )))dt 

K f >() 

+ Cq £ r < X f / (MX^ N + *»V) - MX^ rN ) - a N yd x MX t N ^ N )) ^ 

K ">o Jo J {\v\<r N } v ' \y\ 

+ E E ™L *u (<T N ) - »L-^c m (*&i r N ) 

K N >0 large jump 
1 at C m 

+ E E («L-<„->c» (3-6) 

k n >0 Cm of the 

1 form khjy 

+ E E Km-^L-J^^C?'"") 

K™>0 Cm of the 
1 form (k+l/2)h N 

+ M N - 

Here, the third, fourth and fifth terms correspond to the second term in the right hand side of H3.5|) . 
and Mn is a martingale term given by 



M N = 



.H~„ J "'{|yl<rjv} v 7 



{\y\<rN} 

Equation (|3.6|) can be rewritten 

XjV = + + T^r + T^r + M N , 

where = - , <p ), 1% is the sum of the two first terms in the right-hand-side of (|3.6[) . T% is 
the third one, T$f the fourth one and T% the fifth one. 

The four following Lemmas, whose proofs are postponed to Section T3, 2 1 deal with the asymptotic 
behavior of the terms M N , T% — , and T%. 

Lemma 3.8. It holds 

E\M N \' < 



k'rr 2 r 2 ~ a 



N 
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for some positive constant K. The equivalent term in the case a = 2, 

roc 

M N = cr at y~] / w l Tt d x ip(X^ a )dLl, 
w . „Jo 



satifies the same estimate : 

Lemma 3.9. • It holds 
E 



K?>0 



E\M N \ 2 < 



-T N +T N + r [ L{H*^' rN )d t g t + 2 X ^{H*Vt' rN )d*gt) dt+ [ g Q r,(H * tf> rN )dx 
Jo Ji v ' Jm 

+2c a [ Xt [ [ V(H * jitf ,rjv (x)) (g t (x + a N y) - g t (x) - a N yd x gt{x)) _ 



{iyi<>-jv} w 



iV— >00 



• If r N< \ joN, then 

<Ka%. 



2c a / Xt [ [ V(H * fif' rN (x)) (g t (x + a N y) - g t (x) - o- N yd x g t {x)) 



Uv\<r N } 

The following lemma gives two estimates for the term T N , the first being useful for a constant 
viscosity crjv = o, and the second for vanishing viscosity on — > 0. 

Lemma 3.10. • The error term 



E 



T% + 2c 



\l-X?) [ f v\H*^ rN (x))(H*^(x + o N y)-H*^(x))g t (x)?^ 
JR J{\y\>r N ) V ' \v\ 



vanishes if N 1 r N a goes to 0. 
• It holds 

E\T N \ < K(a NrN - a +a N ). 



Lemma 3.11. One has E\T N \ -> 0. 

N— >oo 



We now have to control the probability for the last remaining term T N to be negative. If 
there is no crossing of particles with opposite signs between fc/ijv and (k + l/2)ftjv, for any k, 
then T N > 0. Indeed, let ^^+V/2)fc 2V < . . . < -X'^+V/^fcjf ^ e a max i m al sequence of consecutive 
particles with same sign. The sequence ^ip(fe+i/2)h N (^(I+i/2)h N )j ls thus a nondecreasing 

sequence, and from the convexity of n and the fact that no particles with opposite signs cross, 
( w (k+i/2)h N )(=!,■■-, 9 i s the nondecreasing reordering of {w^ hN );=i,..., g . Thus, from Lemma 13.131 
below, E K «> fehjv (^(fe+i/2)h„ " ™fch„) ( ^(fc+i/2)h JV (^+ 1 % hjv ) is nonnegative. It is thus sufficient 
to control the probability that two particles with opposite signs cross between fc/ijv and (fc+ l/2)/ijv. 
Since after the murder there is no couple of particles with opposite signs separated by a smaller 
distance than ejv, this does not happen as soon as no particle drift by more than ejv/4 and no 
particle is moved by more than ejv/4 by the small jumps. The drift on half a time step is smaller 
than sup[_j ji \ A'\hN which is assumed to be smaller than ejv/4. We control the contribution of the 
small jumps in the following lemma: 
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Lemma 3.12. Let Bm be the event 

B N = {Vfc < T/h N , Vi, on \^ 1/2)hN - KhZ I < , 

so that no crossing of particles with opposite signs between fc/ijv and (k + l/2)/ijv occurs on -Bjv- 
One /ias, for a < 2, 

/ r^t -a /., \NT/h N 

¥(B N ) > (l-e KhNr x - £ «/ 4CT «''«j 
For a = 2, wie define the event Bn by 

B N = | Vfc < T/h N ,\/i,a N ^L\ k+1)hN - L| fe J < ejv/4} . 

P(Bjv) > (l - Ke" s « /(32hN ' y « ) y T/hN . 
The proof will be given in Section 13.21 

We now gather all the previous information to prove that, depending on the considered case, the 
entropic formulation or the weak formulation holds almost surely. 

1. Constant viscosity <tjv = o , with index < a < 1. 
Define, for Q G V{K), 

F N (Q)= f T)(H*Q )g + [°° f (r,(H * Q t )d t g + 2 X ?i>{H * Q t )d x g) dt 
Jm Jo J& v ' 

+ 2c a [°°(1- X ?) [ [ v'(H*Q t (x))(H*Q t (x + cr N y)-H*Q t (x))g t (x)^^ 
Jo Jw J {\y\>r} \V\ 

+ 2cc f xt [ [ V{H * Qt{x))(g t {x + a N y) ~ g t (x) - a N yd x gt(x))^^^- 
Jo JnJ{\y\<r} \y\ + 

and 

F r (Q) = [ V(H* Qo)go + [ [ (v(H * Qt)d t g + ip(H * Q t )d x g) dt 

\y\l+C 

ixd 

JRJ{\y\<r} " ' 

Notice that from the convexity of rj, one has 

n'{H * Qt(x))(H * Q t (x + ay) - H * Q t (x)) < rj{H * Q t (x + ay)) - r,{H * Qt(x)), 

so that for any < r < r' , it holds F r < F r ' and < F$ . 
From Equation p.6[) . it holds, for N large enough so that rjv < r, 

F N (^ rN ) > F r N N (v N ' r ») = T N + (-Tj, + T 2 N + T N +T A N + M N + F r N N (^)^) . 

From the assumptions made on ejv and Hn one can construct a sequence rjv such that jV _1/ ' a = 
o(rjv), h N rJ, a = o (eivr^ 1 ) and -^e~ eN /iaTN -> 0. Indeed, set r N = e N N~ 1/2X . Then it holds 
N -i/* r -i < KN~ 1/2X and ^r 1 - 01 = h N e- a N {c "- 1)/2X , which vanishes for any value of a. 



-<. j i I V '(H*Q t (x))(H*Q t (x + ay)-H*Q t (x))g t (x)^^ 

JU J{\v\>r] 

+ c a f [ [ ri(H*Qt{x))(gt{x + ay)~gt{x)-ayd x gt(x)) dl '' l '' U 
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Then goes to infinity at the rate of a power of JV, and en/tn = N 1 ^ 2X as well. Thus, 

N e -e N /4crr N vanisneg 
h N 

As a consequence, from Lemmas lT^l3l)ll3T0l and l3TTl E \-T N + T N + T N + Mn + F^ N (fJ, N ' rf 

vanishes as JV tends to infinity, and the event Bn defined in Lemma 13.121 is such that 
P(J3jv) — > 1. On the event Bn, T% is almost-surely nonnegative, so that, from the uni- 
form boundedness of F N with respect to JV, E nN ' rN (F N (Q)~) = E(F N (fJ. N - rN )~) goes to 0. 
To show that the entropic formulation holds almost surely, we need a continuous approxi- 
mation of F N and F r . We define F r,s and Fjj s by replacing every occurrence of H * Qt in 
the definitions of F r and F N by f K l K>t H(. - /(i))7 4 (/(0))d<2(/, k), where 7* is a Lipschitz 
continuous approximation of 7, with P( 7 (Aq) 7^ 7 <s (JVq)) < 5 (see [9], Lemma 2.5, for the 
construction of 7*). Then, for any fixed 8 and r, the family {F r,s } U {F^ s , JV G N} is equicon- 
tinuous for the topology of weak convergence. Indeed, let Q k be a sequence of probability 
measures on K, converging to Q as k goes to infinity. From the continuity of the application 
K, — > K, (/, ft) h-> l K >o/(0), Qo converges weakly to Qo (where Qo and Qo are defined as 
in Q3.2| l1. and from the continuity of the applications K, — > R, (/, k) l K >t7*(/(0))l/(t)<j / on 
the set {(/, k) £ A3, /(t) = f (t— ), f(t) ^ j/}, for all t in the complement of the countable set {t £ 
[0, 00), Q({/(t) ?4 /(«-)} U {ac = t}) > 0}, the quantity / K l K>t H(. - f(t))^/ s (f(0))dQ k (f, k) 
converges almost everywhere to L. l K >t H(.— /(t))7 <5 (/(0))dQ(/, «). From Lebesgue's bounded 
convergence theorem, we deduce that 

sup|F^(Q fc )-F^(Q)| + |F^(O fc )-F^(Q)| -> 

yielding equicontinuity for {F r ' s } U {F^ S ,N £ N}. Moreover, since the sequence \? con ~ 
verges *-weakly to 1/2 in the space L°°((0, 00)), F^ s converges pointwise to F r ' s as JV goes to 
infinity. Ascoli's theorem thus implies that F^ s converges uniformly on compact sets to F r ' s . 
From the weak convergence of iv N ' rN to 7r°°, one thus deduces 

W N,rN [F$ s (Q)-] -> E^[F r ' s (Q)-]. 

N— >oo 

Moreover, for any t > 0, any y, and any probability measure Q satisfying Qo = \uo\ (with Qo 
defined as in (|3.2|) '). which holds true for 7r°°— almost all Q from Lemma 13.51 it holds 

H*Q t (y)- f l K>t #(y-/(*)) 7 4 (/(0))dQ(/, K ) < / | 7 - 7 4 |d| U „| < S, 
Jk Jr 

yielding convergence to for E 7r °°\F r (Q)- - F r - s (Q)-\ +W N ' TN \F N (Q)~ - F^ S (Q)~\ as S 
goes to 0, uniformly in JV. As a consequence, writing 

E^(F r (Qy) <K"°°\F r (Q)- -F r ' s (Q)~\ + \e^ (F r ' S (QD - E""'"" (F^(Q)~)\ 

+ E* N - rN \F r /(Q)- - F N (Q)~\ + E" N - rN (F N (Q)~) 

we deduce that F r (Q) is nonnegative for ir°°— almost all Q. We just have to notice that 
Lemma 13.51 yields that, n°°— almost surely, H * Qo = vo to conclude that the entropy formu- 
lation holds ir°° — almost surely. 

2. Vanishing viscosity on — > 0. 
We define 

Fn(Q)= I v(H*Q )g + I" [ (rj(H*Q t )d t g + 2 X ^(H*Q t )d x g)dt 
Jr Jo Jr v ' 



1G 



and 

F(Q) = / r)(H * Q )go + I I (r)(H* Q t )d t g + i>(H * Qt)d x g ) dt. 



Regularized versions F r / and F s of F^r and F are also considered using the function 7 
instead of 7. In the case a < 2, the same arguments as above, using the second parts of 
Lemmas 13.91 and 13.101 will show that the entropy formulation holds 7r°°— almost surely for 

a 2 r 2 ~ a 

H * Q t , provided there exists a sequence rjv such that — ^ — and a^r]^" vanish, rjv < °n 1 ) 
h N r~ a = o(e iv (<7 J vr J v)- 1 ) and *L e - e »' 4 "' r *' -> 0. 

• For a < 1, any sequence tm vanishing at a very quick rate will fit. 

• For a > 1, since we assumed <tjv < eJ, ° A r ~ 1 ' /A these conditions are satisfied by the 



sequence tm = -^-N 2A (°-i) . 
In the case a — 2, Ito's formula writes 



/ - \ / . \ f( k + 1 ) h N 

<P(k+i)h N [X (k 'l 1)hN j - >fi khN [X kl ;lj = / d t <pt(X^ l )dt 



(k+l)h N 
kh j\ 

f(fe + l)h N . 

+ 2 ! d xVt (X^)A' (H * tf hN (Xf4)) dt 



' kh jy 

2 



+ (Tn I / dl<p t (Xt' l )<lt 

l(kh N ,(k + l)h N ) J{\y\<r} 



+ vn / d x MX"'ldLl. 

J (kh N ,(k + l)h N ) 

The three first terms are treated as in the case a < 2, and the stochastic integral is dealt 
with using Lemma [3.8l For the entropic inequality to holds, we need to control the crossing of 



particles with opposite sign. From Lemma [3.121 if ^-e 32a N h N goes to zero, then no crossing 

occurs. Since our assumptions yield fiNO~% < ej^N -1 ^ and N/hn < KN 1+X , this condition 
holds true. 



3. Constant viscosity crjv = <r, with index 1 < a < 2. 

In this case, since we want to derive a weak formulation, we do not need to consider separately 
large and small jumps. As a consequence it is enough to study the process X^' % . 
Let g be a smooth function with compact support, and define for Q £ V(JC), 



p poo p poo p poo p 

F(Q) = / H*Q o9o + / / H*Q t d t g t dt-a a / / H*Q t (-A)* g t dt+ / / A(H*Q t )d x g t . 
Jr Jo Jr Jo Jr Jo Jr 

Let tp t (x) = f* x gt(y)dy. One has 
1 1 °° 

^ rtXo'^MXo'*) =- x -y( x o' l )f(k+i)h N (x^i 1)hN ) 

,tf>0 fc=0 re iV = ( fe+ i) hjv 

1 00 

+ ^( X o'l (v(k+i)h N (X]l'l 1)hN ) - ip khN (X^ l N )^ 

fc = K, N >kh N 
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From Ito's formula, in the case a < 2, when k\ > kfiN, 



k'+l)h N ) <Pkh N (X kl l ) 



(3.7) 



/ d m (x^)dt + / a^ t (xf- l )A' * ^ hN (x^)) at 

<* [ [ (<Pt{Xt' % + ay) - (p t (X?' % ) - l {ly \< r} ayd x (pt(X" , ' t )) d ^ 

J (kh N ,(k+l)h N ) JI v ' \V\ 

(MX^ + ay) - MX?-*)) *(dv,dt). 



+ c, 
+ 



(kh N ,(k+l)h N ) 



We denote r t = max{fc/ijv, fc/ijv < t}. Multiplying (|3.7|l by -fil K N >kh f(X summing over : 
and k, and integrating by parts, one obtains 



IS ~" 



9t3«H * At d* - 



(-A)* flt tf*/if dt 



+ 



£ 7«'*) (pttC'* + ^) - ^ t (Xfj 1 )) u\dy,dt) 



(0,oo)xK „ 



1 °° 

+ l)hiv(^(fc + l)h jv )' 



TV 



(3.8) 



fe=0 K? = (k+l)h f 



Combining an adaptation of Lemma 13.141 stated in Section 13.21 with A replacing 77, and 
integrating by parts, the difference 

jr r E l(X?' i )d x MX?' i )A'(H*il khN (X^l))at+ r f d x9t A(H*^)dt 

vanishes in L 1 . Using an adaptation Lemma 13.81 the the fourth term in the right hand side 
of (13.81) vanishes in L 2 . The fifth term vanishes in L 1 since 



1 00 

jvE E 7(^f ,l )¥'(fc+i)h J v(^(Iii)h J , 



^f=(/ ! +l)^j v 



-jvE E |^(*-'+i)hjv (-^(I+ijhjv) - <P(k+i)h N {x^i 1)ht 



fc— pairs killlcd 
at time (k-\-l)hjy 

^Ken- 
As a consequence, E" |_F(Q)| = E\F(fi N )\ vanishes. We conclude by regularizing the func- 
tion 7 as in the two first points, that \F(Q)\ = 0. Thus, F(Q) = almost surely, so that 
H *Q almost surely satisfies the weak formulation. 

The case a — 2 is treated in the same way, the only difference lying in the nature of the 
stochastic integral. 
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3.2 Proofs of Lemmas Q to I3TT21 

In this section, we give the proofs of the previously admitted lemmas of Section [3T 



Proof of Lemma \y.8\ Since the particles are driven by independent stable processes and since the 
inequality \wl\ < It holds for some constant K not depending on t, i and N, 



.n >0 J o \J{\y\<r N } v ' 

<2a 2 NCa E [ V /°°«) 2 xf / 
W>° ° J l 



I I, H n2 dydt 



fll/l<^JV> 



■ a N r N 



<K~-^-] o listed*. 

A similar proof with stochastic integrals against Brownian motion yields the result for a = 2. □ 
Proof of Lernma \3.9\ Integrating by parts, one finds 

N roc . roc r roc p 

£ / dt = - / n(H* H?' TN )d t g t dt + / / n(ji?> r »(m))d t g t dt 

i=1 JO JO J8 JO Jl 

yielding, from Lemma 13.141 below, 

£ r<d m (x t N —)dt+ r [ r,(H*il?> r »)d t g t dt- [°° [ v (^> r 
i=1 Jo v ' Jo Jr Jo Jr 



gtdt 



— > 0. 

N— J-oo 



From the constancy of /i^ ,rjv (lR) and an integration by parts, one has 

—Tn + / / v(^' rN mdtgtdt = - / g 0V (H * 
Jo Jr Jr 

Another integration by parts yields 



N 

2c, 



~[Jo J{\y\<r N } v / \ / \ 

/ Xt [ [ (gt(x + a N y) -g t (x) - a N yd x g t (x)) r/(H * ^' rN (x))^^- 

Jo J{\y\<r N } Jr 13/1 

: Xt r l{i l 't' rN W) / / (fft(* + onv) - - ^ya. ft (x)) 

Jo J{|»I<i-n} Jk ^ 

/ Xt f [ (gt(x + a N y) -g t (x) - a N yd x g t {x)) r/(H * fj,^ rN (x))^^ 

Jo J {\y\<r N } Jr \V\ 



= -2c, 
+ 2c, 

= - 2c, 

Moreover, from the regularity of A and 77, it holds 



dydt 

dxdydt 



dxdydt 

+a ' 
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so that 

JV 



E 



from an adaptation of Lemma 13.141 (replacing 77 by ?/> in the definition of w l t ). This concludes the 
proof of the first item of Lemma [3.91 

To prove the second item, observe that the change of variable z = umV yields, for rjv < — , 



2c c 



<2c a (T% 



Xi 



{\v\< t n} 



T)(H * p,?' rN (x)) (g t (x + a N y) - g t (x) - a N yd x g t (x)) "' 



N 

Xt 



\v(H * h ,r,v (^)) (9t(x + z) - g t (x) - zd x gt(x)) 



dzdxdi 



{|*l<i} 



□ 



Proof of Lemma \3.1(A First notice that 

E fVxf) / 



'>0 



{|vl>nv} 



with p defined by the following formula: (/ij /,J ' iv ' rjv being the measure obtained by moving in the 
expression of fi t ' rN the particle X^''' ,N to the position X t ' 8,rw + omv) 



To prove the second item in Lemma 13.101 we integrate by parts, and, using the definition 
of p l v ''' N ' rN and the compactness of the support of g, it holds 



/ Vt A P r = / g t ( V (H*fl> i ' N ' rN )-v(H*ii?' rN )) 
Jm Jm v ' 



< K 



(o-Ny) A 1 
N 



(3.9) 



so that 



nn\<K (l-x* 



dydt 



) {onv) A 1 < K(a% + a N r N °). 

J{|a/|>rAr} 12/1 



Now let us prove the first item of Lemma 13.101 Applying the same martingale argument as the 
one used to prove E|M;v| — ¥ 0, and using the upper bound K/N in (|3.9p . one has 



T A N -2c a (1-xf) 



■/{i»i>r„} v.ft. ' j \y\ + 



K 



Let us give a more explicit expression for p\ '*. For simplicity, we denote 



I 



It E 7(xS)i 



/ 
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and for i 7^ j, 



*i J ' ± = 1 ««>t 1 ««>t 



One can write 



h E 7(^)1, 



AT 



AT 



* E 7(^)1, 



N 



Mi 



E ^ = E •"'^f-' 

+ E E 
/ 



E ^ 



•toi 1 



+ E E(^' 



<<X™- ] - r "+<r N y 



(3.10) 



In this expression, the two first terms deal with particles jumping from the site X t ' l ' TN to the 
site X t *' rN + &Ny, while the third term corresponds to the jump from right to left of the particle 
labelled j above the particle labelled i and, conversely, the fourth term corresponds to the jumps of 
particle j from left to right over particle i. Notice that this last equality, as well as (|3.11[) below, 
only holds when each X^ ,l ' TN + o-jvj/ is distinct from all X^^' rN . However, for all t, this condition 
holds dy-almost everywhere, which is enough for our purpose. 

In the entropic formulation (|2.3[) . the term that should appear for large jumps is given by 



2c D 



J{\y\>r N } 



iptdal 



dydt 



where 



r? =d x (r,'(H * Af") (H * ^- rN (- ~ °NV) - H * ^' rN )) 



4 E 7(^)'/'(ff*Af ,rjv (^ Ar,< ' rw +^y) 

K i>* K i>* 

+ Y, ( H * Af" - o N y) - H * Af' rw )) 

K *> >t 

x (r, 1 (H * (X^ r » )) -r,'(H* tf' r » (Xf'^ -)) ) 5 xf 



I Y 7(XiW(H*^(X 



(3.11) 



When computing the difference p\ — a\ integrated against some bounded function, using Taylor 
expansions for r), one can check that, up to an error term of order O (jr) the first terms in the right 
hand side of (|3.10|l and (|3.11|l cancel each other, the second terms as well, and so does the sum of 
the two last term in (|3.10p with the last one in (|3. 1 1 [) . Consequently, 



r< 1 - *n 1 if'***) - r (i - ^ 1 if***) 1^ 

Jo J{\v\>r N } WR / \V\ JO J{\y\>r N } \J K / \V\ 



K 



This concludes the proof. 



□ 



21 



Proof of Lemma \3.11\ For a time £ m of the form fchjv, no particle moved in the interval (Cm-i, Cm), 
so that w\ m — _ — 0, unless the particle labelled i has been killed at time £ m . Hence, 

JV 

T ^=E E ( w L-™L->< m (xt i,rN ) 

i=l Cm of the 
form khjy 

=- E E »L-^c m «r N ). 

Cm of the K N = Q m 
form fc/ijY 1 

This sum is actually a sum over pairs of close particles with opposite signs, thus 



Cm of the pairs of particle 

form fc/ijv killed at time Ctt 



E 



Cm. of the pairs of particles 

form feft.^ killed at time Cm 



Indeed, a couple of killed particles is such that \X^' rN — X^' rN \ < £jv and is made of 

particles with opposite signs, so that 



( 7 (AS) + tW)) r,'(H * (Xf ' ^)) + O (±^j 



K 
~N2- 



□ 



Proof of Lemma \3. IB Notice that from independence of the increments, denoting by £- r a Levy 
process with Levy measure c a l^ y \ <r . fey , it holds 



at) =P(<t^7| <e N H) NT ' h » 



= 1-P ajvnvl^ _ a | >£iv/4 



NT/h N 



Since the Levy measure Cal|y|<i i 8 a+a nas compact support, the random variables Cf 1 have expo- 
nential moments, and Chernov's inequality yields 



a N r N \Cf L _ a \ > ejv/4 < E e 

h Nr N 



h N r N " I g- e w/ 4o "WN _ e A 'hW Jv a - e w/ 4o "WN 



where the constant K does not depend on N. 

In the Brownian case a = 2, we use the tail estimate e -1 da; < Ke~ M for positive M. □ 

Lemma 3.13. £ef a\ < ... < ajv and bi < . . . < 6jv be two nondecreasing sequences of reals 
numbers. Then the quantity aib a ^ for some permutation a is maximal when o(i) = i for 

all i. 
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Proof. From optimal transportation theory (see [141 page 75]), the quantity X^ili( a ! — ^o-(i)) 2 * s 
minimal when a is the identity. Expanding the square, we see that X]^Li( a » — °<r(i)) 2 = X^ili( a ? + 
6?) - 2 2i=i a ibcr{i)- Thus, ^i=i a i & o-(i) is maximal if and only if J2i=i( a i - °<t(i)) 2 is minimal, 
concluding the proof. □ 

Lemma 3.14. Let f be some bounded function with compact support on [0, oo) x R which is smooth 
with respect to the space variable. If /ijv vanishes and o~n is bounded, it holds 



lim E 

N— >oo 



e r(w-<)/ t (xf^ 



'>0 



(If 



Proof. First notice that when t is not in an interval [fc/ijv, (k + l/2)/ijv], it holds = w\ t , since no 
particle moved between n and f. Then, one can write, from the assumptions on /, 



Kf >0 



Integrating by parts, it holds 



' >o 



■»>o Jo V ' 
^ xf f (v(H* Vt' rN (*)) - i? (# * An' rN (*)) ) ft./t (^)d^df 



/ N 



E 1 X^ i ^<x<xi 
K N >0 



+ 1 x"- i -'«<x<x t JV ' l - r « J ^/tO^dadt 



^§ fx? E kf^-^^lAidt 

K f>o 

We conclude the proof by writing 

E ^ ^ Xf E - | A ld * = K jf X f 1 K » >0 |*f " 



' N A ldt 



< T I h N sup +E ((ff W |A^ 1 ' r "|) A l) 



This last quantity vanishes when ftjv goes to 0. 



□ 



4 Numerical results 

In this section, we illustrate our convergence results by some numerical simulations. We simulated 
the solution to the fractional and the inviscid Burgers equations 



dtu+-d x (u) + a a (-A)- =0 and d t u + -d x (u) = 0, 



23 



corresponding to the choice A(x) = x 2 /2, with different values for the parameter a. 

One can find an explicit exact solution to the inviscid Burgers equation (see [T5]) and we compare 
the result of the simulation to this exact solution in the vanishing viscosity setting. However, to our 
knowledge, no explicit solutions exist in the case of a positive viscosity coefficient for a < 2, so that 
we have to compare the result of our simulation with the one given by another numerical method. 
Here, we use a deterministic method, introduced by Droniou in [7]. 

4.1 Constant viscosity (a^ = cr) 

We give three examples of approximation to the viscous conservation law. On Figures [T] [5] and 
we show the approximation of the viscous conservation law with respective index a = 1.5, a = 1 
and a = 0.1 and diffusion coefficient a — 1 using N = 1000 particles, with parameters h — 0.01 and 
£ = 0.04 at simulation times 0.25, 0.5, 0.75 and 1. The continuous line is the simulated solution, 
and the dotted line is the "exact" solution obtained with the determistic scheme of [7] using small 
time and space steps. 

We now investigate the vanishing rate of the error, that is the Riemann sum on the discretization 
grid associated to the integral in Theorems l3.ll I3^2l and l3^3l On Figure [4] is depicted the logarithmic 
plot of the error as a function of N where we used the relation /ijv = W/N, and £jv = 40/iV, with N 
ranging from 10 to 10000, in the three cases a = 0.5, 1 and 1.5. In the case a < 1, this relation 
between N, /ijv and £jv satisfies the condition of Theorem 13.11 These pictures make us expect a 
convergence rate of -^=, corresponding to the optimal rate analyzed theoretically in [510, in the 
case a = 2, without killing. 




-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 



Figure 1: Approximation of the conservation law with index a — 1.5. 



4.1.1 Behaviour as h — s> 

We give in Figure [5] the approximation error at fixed number of particle, with a vanishing time 
step h, in logarithmic plot. We set the parameter e to be equal to 4/i so that the condition of 
Theorem [3T] is satined. We took N = 340000 and a = 1. We set a = 0.5, a = 1 and a = 1.5 
respectively. The different parameters h range from 1 to 2 -8 . In it is shown, in case a = 2 and 
the initial condition is monotonic, that the error is of order h. In view of Figure [5j it seems that the 
convergence rate is still of order h, even for a < 2 and any initial condition with bounded variation. 
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Figure 2: Approximation of the conservation law with index a = 1. 




Figure 3: Approximation of the conservation law with index a — 0.1. 

4.2 Vanishing viscosity (a N — > 0) 

We consider the Burgers equation 

d t v = d x (u 2 /2) 

with initial condition uo(x) = lr_ 3i _ 2 ] — 1[2,3], which is the cumulative distribution function of the 
measure <5_3 — <5_2 + <$2 — £3. In that case, the solution of the Burgers equation is explicit and given 
by the expression 

/x + 3\ fx — 3 \ 

u(t,x) =min y— — ,1 1 l [ _ 3 , min( _ 2+ t,_ 3+ V2t,o)] +max( — — ,-1 J l [roa x(2-|,3-^,o),3]- 

We compare the function u to the function obtained by running the Euler scheme with a small 
diffusion coefficient a. One can expect the approximation to be better for large values of a. Indeed, 
for small values of a, the particles tend to jump very far away, and subsequently "disappear" from 
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Figure 4: Logarithmic error in the approximation of the conservation law with index a — 0.5, 1 and 1.5. The 
respective slopes are —0.46, —0.41 and —0.56. 




-5 -4 -3 -2 -1 -0 -6 -5 -4 -3 -2 -1 -0 




-6 -5 -4 -3 -2 -1 -0 



Figure 5: Logarithmic plot of the error as h vanishes, with a fixed number of particles, at respectively a = 0.5, 
1 and 1.5. The slopes are equal to 1 up to an error of 0.01. 



the simulation. The consequence of this behavious is that the solution is somehow decreased by a 
multiplicative constant. 

For large values of a, the approximation is quite good, even for not so small diffusion coefficients. 
Figure [6] gives the result of the simulation of the Euler scheme with parameters a = 1.5, e = 0.04, 
a = 0.1 and h = 0.01, at the different times 2, 4, 6 and 8 for N = 10000 particles. Figure [7] gives 
the same simulation for a = 1. In the case a < 1, and especially when a is small, one need to take 
a very small value for the diffusion coefficient in order to have a reasonable approximation of the 
solution. Indeed, the approximation depicted on the Figure [S] is the appproximation of the solution 
at times 2, 4, 6 and 8 for diffusion coefficient a — 10 -4 . Here, we used 10000 particles killed at a 
distance e = 0.01, the time step being h — 0.01. On Figure [9] we show the same simulation, with 
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diffusion coefficient changed to a = 10 




Figure 6: Approximation of the inviscid conservation law by a fractional Euler scheme with index a = 1.5 and 
diffusion coefficient 0.1. 




Figure 7: Approximation of the inviscid conservation law by a fractional Euler scheme with index a = 1 and 
diffusion coefficient 0.1. 
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